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ABSTRACT i 

Pool boiling heat- transfer measurements were made using 
seven Gewa-T copper tubes in R-113. The first tube (19-mm 
OD xQ.75 fins/mm X 0.25 mm gap) was tested under three condi- 
tions: (a) plain, (b) with four shrouds and (c) with up to 

5 wires (0.13 mm diameter) wrapped in the inter-fin trough 
area. Verifying prior data, the shroud with longitudinal 
openings of 60° at the top and 8.5° at the bottom gave the 
best performance among the four shrouds tested. This shroud 
increased the boiling heat-transfer coefficient by 253 per- 
cent (over the smooth-tube case) at a heat flux of 10,000 
2 

W/m while it showed only a 18 percent increase at 100,000 
2 

W/m . When wire wraps were provided, in all cases, the 
heat- transfer coefficient was improved at all heat fluxes. 

The use of three wires gave the best performance with 341- 

percent increase in heat- transfer coefficient at 10,000 

2 2 
W/m and a 130 percent increase at 100,000 W/m . 

The next three Gewa-T tubes had dimensions of 25 mm 

OD xO.75 fins/mm with fin-tip gaps of 0.15, 0.25 and 0.35 mm, 

respectively. The last three Gewa-T tubes had the same 

dimensions except that the fin density was 1.02 fins/mm. 

For each of these fin densities, the 0.25 mm fin-tip gap 

produced the best boiling performance at all heat fluxes. 

Also, for a given fin-tip gap, the boiling performance 
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increased with increase in fin density. For a fin-tip gap 
of 0.25 mm, the lower fin-density tube produced an 80-percent 
increase, while the higher fin density tube produced a 103- 
percent increase in the boiling heat-transf er coefficient. 
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INTRODUCTION 



I . 



A. BACKGROUND 

The energy crisis that began in 1974 has made a consider- 
able impact on the design and operation of energy-conversion 
devices. Today, designers are challenged to design more 
and more efficient heat-transfer systems. Among the numer- 
ous areas that designers are faced with this challenge, 
boilers or evaporators used in many engineering applications 
demand considerable attention. The use of advanced boiling 
surfaces may lead to significant improvements in efficiency 
and/or significant savings in the initial cost and equipment 
sizes. Also, due to the ever-decreasing sizes of electronic 
components, the presence of very high heat fluxes has made 
the cooling of such components a challenging problem. Marto 
and Lepere [Ref. 1] stated that the use of advanced boiling 
surfaces may have potential promise in the field of electronic 
cooling. However, the lack of generalized performance data 
and useful theoretical models for such advanced surfaces 
have still kept designers away from the utmost use of these 
surfaces . 

Due to a lack of complete understanding, the theoretical 
treatment of the nucleate pool-boiling process, even on a 
smooth surface , is extremely complex. The existence of a 
large niimber of semi-empirical formulas is probably a very 
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good indication of the complex nature of the nucleate boil- 
ing process. Chongrungreon and Sauer [Ref. 2] summarized 
nine(9) correlations for nucleate pool boiling of single- 
component fluids from smooth tubes. They stated the inade- 
quacy of any of these correlations when compared with 
experimental data. 

Nucleate pool boiling is strongly dependent on the sur- 
face structure; i.e., the heat flux increases with increase 
in stable nucleation site density. The size, shape and den- 
sity of such nucleation sites are almost impossible to pre- 
dict. However, the entire nucleate pool-boiling process is 
dependent on the performance at these nucleation sites. It 
is we 11 k nown that the ability of these sites to develop and 
trap small vapor bubbles is the reason for improved performance 
in nucleate pool boiling in comparison with natural-convective 
boiling . 

When a bubble is being generated at a surface nucleation 
site (or a pore) , the wall superheat is given by the follow- 
ing approximate theoretical equation [Ref. 3]: 

2 0 T„ 

T - T = — 

W S hpg R 

Thus, for a given fluid, the wall superheat can be de- 
creased (boiler performance increases with decrease in wall 
superheat) with increase in pore radius. Of course, this 
radius must have a practical upper limit since too large a 
pore would not trap a vapor bubble. 
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Nishikawa and Ito [Ref. 4] recently discussed two methods 
to augment nucleate boiling: (1) use teflon-coated pits to 

reduce wettability of boiling liqud, and (2) manufacture 
surfaces with numerous re-entrant cavities that have the 
ability to trap vapor and keep the nucleation sites active. 
Due to the high wettability of low surface-tension fluids 
(such as f luorcarbons) , the use of teflon-coated pits has not 
been successful. Thus, the surfaces with re-entrant-type 
cavities are the viable candidates for low surface-tension 
liquids . 

Taking advantage of the concept of re-entrant cavities, 
a number of enhanced boiling surfaces have been invented 
during the last two decades. They include the Linde High 
Flux surface of Union Carbide, Gewa-T surface of Wieland 
and Thermoexcel-E surface of Hitachi. All of these surfaces 
are provided with re-entrant cavities through specialized 
manufacturing processes. These re-entrant cavities introduce 
yet more parameters to be desired. 

In view of the inadequate understanding of the nucleate 
pool boiling process on smooth tubes, a thorough understand- 
ing of the boiling process on advanced surfaces may take 
many more years to come. Thus, pbesent-day designers have 
very little information when considering advanced boiling 
surfaces. It is very important that researchers gather a 
vast amount of experimental data covering a wide range of 
heat-transfer fluids. Such extensive data will enable the 
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development of generalized correlations or successful 
theoretical models. 

Yilmaz, et al. [Ref. 5] compared the nucleate pool-boiling 
heat-transf er performance of three copper tubes: Linde High 

Flux, Wieland Gewa-T and Hitachi Thermoexcel-E . Their data 
taken with p-xylene showed that while all tubes performed 
better than smooth tubes, the best performance was with the 
Linde High Flux surface. The Thermoexcel-E tube performed 
better than the Gewa-T tube at low heat flux and both per- 
formed similarly at high heat flux. These data were taken 
under the fully-established boiling regime; thus, no infor- 
mation was provided on the hysteresis behavior or temperature 
overshoot of these surfaces. 

Bergles and Chyu [Ref. 6] compared the nucleate pool- 
boiling heat-transfer performance of four different Linde 
High Flux tubes to a plain copper tube in distilled water and 
in R-113. They used three different methods of aging the test 
surface prior to collecting data. They reported significant 
temperature overhsoots of both plain and High Flux surfaces 
in R-113. They also reported that these overshoots were 
sensitive to aging, initial subcooling and rate of power 
increase . 

Marto and Lepere [Ref. 1] obtained results on the heat- 
transfer performance of three heat-transfer surfaces: Linde 

High Flux, Thermoexcel-E, Gewa-T, and a plain copper sur- 
face in R-113 and FC-72. They reported a two to tenfold 
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increase in the heat-transfer coefficient when compared to 

plain surface. The High Flux surface was most effective over 

a broad range of heat fluxes. The Thermoexcel-E surface 

showed similar gains in heat-transf er coefficient to that 

2 

of the High Flux surface below 10 kW/m . The Gewa-T surface 

was not as effective as the othe'r surfaces at low heat fluxes 

2 

but performed comparably at high heat fluxes (near 100 kW/m ) 
Marto and Hernandez [Ref. 7] studied the boiling per- 
formance of a Gewa-T finned surface in Freon-113. They 
examined the influence of placing aluminum shrouds around the 
test tubes, and found that the thermal performance of these 
surfaces were dependent upon the liquid-vapor flow within 
the channels between neighboring T-shaped fins. They also 
examined the influence of the fins by progressively machining 
away the fin height to arrive at a smooth cylindrical surface 
They reported that the addition of the "T-caps" to straight 
fins produced the most significant improvement in heat- 
transfer performance when compared to a smooth tube. The 
use of shrouds increased the heat- transfer coefficient at 
all practical heat fluxes, but decreased the heat- transfer 
coefficient at very high heat fluxes when compared to the 
unshrouded surface. 

B. THESIS OBJECTIVES 

Based on the foregoing discussion, the major motivation 
in this work was to further study the boiling heat-transf er 
performance of Gewa-T surfaces in R-113, and to attempt 
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to improve this performance through the use of further 
enhancing techniques. 

There were four major objectives in this thesis study: 

1. Take data to verify results of Hernandez [Ref. 8] on 

a standard 21.2 mm O.D. Gewa-T tube with and without shrouds. 

2. Use wire wraps in- the inter-fin area to improve 
nucleate-boiling heat-transf er performance from Gewa-T 
tubes . 

3. Take data using the combination of best shroud (found 
in 1) and the best number of wire wraps (found in 2) on 
Gewa-T tubes. 

4. Take data to study the effect of fin-tip spacing and 
fin density on the heat- transfer performance of Gewa-T 
tubes . 
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II. EXPERIMENTAL APPARATUS 



A. REQUIREMENTS GOVERNING DESIGN 

The purpose of this research was the observation of 
nucleate pool boiling from different Gewa-T surfaces. In 
order to get a more comprehensive knowledge about the physi- 
cal mechanisms that govern this phenomenon, it is necessary 
to measure the parameters involved which affect the nucleate 
boiling process. The selected parameters to be obtained 
are the following; 

1. Barometric pressure 

2. Test section wall temperature 

3. Fluid bulk temperature 

4. Vapor temperature 

5. Surface heat flux 

B. TEST APPARATUS 

The experimental apparatus designed by Lepere [Ref. 9], 
and modified by Hernandez [Ref. 8], with further minor modi- 
fications, was used during this thesis work. To gather data 
more conveniently at a faster rate, a Hewlett Packard model 
9826 computer and an HP-3054A automatic data-acquisition 
system were added to the experimental facility. Lepere 
[Ref. 9] and Hernandez [Ref. 8] provide detailed descriptions 
of the apparatus; therefore, only a brief description is 
presented in this thesis. 
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Figure 1 shows a schematic of the test apparatus. The 
boiling surfaces consist of a cylindrical copper tube with 
a smooth interior surface and a commercially-available T- 
shaped finned surface manufactured by Wieland-Werke , AG. Each 
tube was conditioned for the insertion of thermocouples by 
drilling eight longitudinal holes as indicated in Figure 2 
and introducing the respective copper-constantan thermocouples 
centered by means of capillary copper tubes. In the center 
of the tube a cartridge heater was placed to act as the heat 
source as shown in Figure 3. 

Both ends of the test section were machined to remove 
the enhanced surface on the outside and to increase the 
interior diameter. This reduced the cross sectional area 
and thereby minimized the heat transfer from the ends of the 
test section. A Teflon plug was inserted into the free end 
for sealing and heat transfer reduction purposes. The other 
end of the test section was connected to a long copper ex- 
tension tube through a 90-degree elbow, which allowed the test 
section to be placed in a horizontal position within a thick- 
walled Pyrex glass vessel. Since one of the requirements 
of the apparatus was to maintain the liquid at saturation 
temperature and atmospheric pressure, a plate heater and a 
vent tube were provided to satisfy these conditions. In 
addition avoiding excessive loss of vapor, a water-cooled 
copper coil was placed in the vapor space of the vessel, to 
act as a primary condenser. The vessel was covered with a 
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Figure 1. Schematic of Test Apparatus 
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Figure 2. Thermocouple Radial and Angular Positions 



CEWA-T SURFACE SOLDER THIN WALL 
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Figure 3'. Cross-Sectional View of Test Section 




Plexiglas lid, which had been fitted with a rubber 0-ring. 
Several holes were drilled in the lid to fix in place the 
test section extension, vent tube, connections for primary 
and secondary condensers, and two more thermocouples to sense 
both the bulk liquid and vapor temperatures . The Pyrex 
glass vessel with the attached fittings was placed on the 
plate heater, which served to maintain the liquid at satura- 
tion temperature. A secondary condenser was mounted external 
to the vessel. 

C. INSTRUMENTATION 

The power supply was connected to the cartridge heater 
through a voltage regulator and then to a Variac scaled from 
0 to 120 VAC. This arrangement allowed fine adjustment for 
a stable feed voltage during the runs. A precision resistor 
of 2.031 ohms was connected in series with the cartridge 
heater to provide a measure of the current flow by reading 
the voltage drop across this element. A digital voltmeter 
was used to measure both the voltage across the resistor and 
the voltage across the cartridge heater. All the thermo- 
couples were connected to a computer data acquisition system, 
and ambient conditions and voltage readings were supplied 
manually. The microcomputer processed all this information 
through a program included in Appendix A. 

D. CALIBRATION OF THERMOCOUPLES 

A total of fifteen wall thermocouples in tubes 1 and 7 
(one thermocouple in tube 1 was defective, hence it was 



24 



excluded) , and the two thermocouples for vapor and liquid 
temperature measurements were calibrated in an insulated, 
stainless-steel bath. The tubes were completely immersed 
into the bath, which was provided with a motor-driven mixer. 

The bath temperature was measured using a platinum-resistance 
thermometer accurate to ±0.01 K. 

Since the temperature measurements during this investiga- 
tion ranged from 45°C to 60 °C, thermocouple calibration was 
performed in the range of 45°C to 70°C. To begin the cali- 
bration process, hot water (at about 75°C) was added to the 
insulated container. Following proper mixing, all thermocouple 
readings were automatically recorded through the data- 
acquisition system. The bath temperature was measured using 
the platinum- resistance thermometer. To obtain other lower 
temperatures, small quantities of water were added to the 
bath, allowing 5 minutes for temperature stabilization of the 
test section. 

The discrepancy (i.e., the platinum- resistance thermometer 
reading minus the thermocouple reading) was plotted against 
the thermocouple reading as shown in Figure 4 . The second 
order polynomial curve shown in this figure was generated 
using the data from all seventeen thermocouples. However, 
for clarity, only three thermocouples were selected for data 
points from each tube. It is evident that the thermocouples 
in both tubes show similar discrepancies, the thermocouple 
readings have only about ±0.1 K discrepancy around the 
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Figure 4. Thermocouple Calibration Curve 



calibration line. Based on this fact, no calibration was 
performed for the other five tubes. Instead, the calibration 
curve generated for tubes 1 and 7 was assumed for all tubes. 
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III. EXPERIMENTAL PROCEDURES 



A. NORMAL OPERATION 

The Procedure B employed by Marto and Hernandez [Ref. 7] 
has been demonstrated to be the most adequate to carry out 
the boiling runs. Because of the large influence of surface 
past history on boiling incipience, this procedure was estab- 
lished to allow repeatability and permit comparison of the 
obtained data for different enhancements. 

The procedure is as follows : 

1. Once the test section has been immersed in the pool, 

it was subjected to one hour of pre-boiling by setting the 

2 

cartridge heater to give a heat flux of about 30 kW/m . The 
plate heater was adjusted to the maximum temperature for this 
one-hour period. This initial vigorous boiling served to 
degas the liquid and to force the noncondensable gases out to 
the atmosphere. 

2. After the initial aging process, both test-section heater 
and plate heater were secured, allowing cooling of the liquid 
and the test section for about 30 minutes. 

3. The plate-heater voltage was then adjusted in order to 
maintain the liquid at saturated conditions and the cartridge 
heater was turned on to the first setting of 12 volts. 

4. For all the consecutive settings during both increasing 
and decreasing heat fluxes, the boiling was allowed to 
stabilize for five minutes at each power setting. 
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5. At a given power setting, the following data were 
recorded: heater voltage, precision-resistor voltage, vapor 

temperature, liquid bulk temperature and wall temperatures 
of the test section. 

B . GEWA-T SURFACES 

Three groups of surfaces were tested. The f irst . surf ace 
was the same as used by Hernandez [Ref. 8] with a fin density 
of 0.75 fins/mm, an outside diameter of 21.2 mm, and an 
inter-fin gap of 0.25 mm. The second group of three tubes 
had 0.75 fins/mm, an outside diameter of 23.1 mm and inter- 
fin gaps of 0.15, 0.25 and 0.35, and the third group of three 
tubes had 1.02 fins/mm and all other specifications were the 
same as the second group. 

The geometry of the enhancement of these seven Gewa-T 
surfaces was observed by making a longitudinal cut of the 
tubes and photographing them with 50X magnification (Figure 5) 
As can be seen, the tubes with 0.75 fins/mm had an approxi- 
mately circular cavity shape, while the 1.02 fins/mm tubes 
had a somewhat elongated cavity shape. This difference in 
shape was probably due to the technique used in construction 
to obtain the various specifications. All the T-shaped fins 
had a fin height of approximately 0.95 mm. 

The exact determination of diameters at the base of the 
fins was carried out from the magnified photographs of the 
cavities. This diameter was determined by subtracting twice 
the fin height from the measured outside diameter: 
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(a) Gewa-T Tube Number 2 



(b) Gewa-T Tube Number 5 



Optical Micrograph of Gewa-T Fin Profile, SOX 



Figure 5. 
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O.D. 



D 



2 



- 2Hp , 



resulting in the dimensions presented in Table 1. 



TABLE 1 

Dimensions of Gewa-T Finned Tubes 



TUBE 

NO. 


FINS/ 

mm 


FIN TIP 
GAP (mm) 


FIN HEIGHT 
(mm) 


O.D. 

(mm) 


D2 

(mm) 


1 


0 . 75 


0 . 25 


0 . 865 


21.2 


19.47 


2 


0.75 


0.15 


0.980 


24.99 


23.05 


3 


0.75 


0.25 


0.960 


24.99 


23.07 


4 


0.75 


0 . 35 


0.970 


24.99 


23.05 


5 


1.02 


0.15 


1.0 50 


25.30 


23.20 


6 


1.02 


0 .25 


1.050 


25.30 


23.20 


7 


1.02 


0.35 


1.125 


25.30 


23.05 



For simplicity, the diameter at the base (D2) was assumed 
to be 19.5 mm for Tube 1 and 23.1 mm for all other tubes. 

C. HEAT-FLUX CALCULATION 

According to earlier description about test-section con- 
struction, the cartridge heater acts as the heat source, 
inserted into the tube, and variations in heat-transf er rate 
are made through different power settings with a variac. 

= VI (1) 
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where ; 



Q„ = heat-transfer rate from the cartridge 
heater (W) 

V = voltage across the cartridge heater 

element (volts) 

I = current through the heater element (amps) 

This circuit is connected in series with a precision 
resistor, which provides the means to compute the electrical 
current by the measurement of the voltage drop across the 
resistor : 

= ''res^'^ <2) 

where : 

V = voltage drop across the precision 

resistor (volts) 

R = resistance of the precision resistor (ohms) 

= 2.031 ohms 

Both ends of the Gewa-T tube were machined to provide a 
smooth surface and a reduced wall thickness in order to mini- 
mize longitudinal heat conduction and therefore the heat 
losses from the unenhanced ends. The geometry of these ends 
is shown in Figure 6, where: 

b = base of the T-shaped fins 

T, = temperature at the base of T-shaped fins 
^ (°C) 
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Lu , 



Figure 6 . Geometry of Unenhanced Ends 



For this cylindrical geometry: 



A 



C 




( 3 ) 



A_ = Px L 
S c 



( 4 ) 



where : 

2 

A^ = cross-sectional area of tube wall (m ) 

= tube outside diameter (m) 

D. = tube inside diameter (m) 

1 

2 

Ag = tube outside surface area (m ) 

P = tube outside wall perimeter (m) 

It was assumed that the temperature at the base of the 
thin-walled fin was equal to the temperature at the base of 
the T-shaped fins. Therefore: 
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0(x) = T(x) - 



(5) 



In general. 




2 



m 



6 = 0 



( 6 ) 



where : 



2 



m 



h P 
kA^ 



(7) 



Therefore 



Qf 



Qb 



-kA 



dT 
C dx 



x=0 



-kA. 



dx 



x=0 



( 8 ) 



Now, assuming negligible convective heat loss from the 
thin tube tip: 



de 




0 



(9) 



and 



6 _ Cosh m(L-x) 

6,3 ~ Cosh mL 



For this case 



Q 



F 



vt 



h F kc Ac 



0g Tanh mLj^ 



( 11 ) 
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6 (o) 



( 12 ) 



e 



B 



B 



- T 



SAT 



where : 



Q 



F 



T 

SAT 



k 



c 



h 



heat-transfer rate through unenhanced 
end (W) 

fluid saturation temperature (*^C) 

thermal conductivity of the copper (W/m.K) 

2 

average heat- transfer coefficient (W/m . K) 



Assuming the wall temperature is the same as the average 
of the eight calibrated wall temperature thermocouple 
measurements, less temperature drop due to radial conduction 
from thermocouple position to the base of the T-shaped fin, 



ln(D«/D ) 

T = T - O = — 

B AVG 2 7T L 



(13) 



T 

AVG 



8 

( I T )/8.0 
n=l ^ 



(14) 



where : 

= diameter at the position of thermocouples 
(m) 

Tavg “ average wall temperature (°C) 

For the geometry shown in Figure 6 : 



L 



C 




(15) 
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t 



D. 

1 



( 16 ) 



= D - 



where : 

= corrected length (m) 

Ly = length of unenhanced ends (m) 
t = tube thickness (m) 

The average difference between the wall temperature and 
the saturation temperature may now be determined from the 
following equation: 
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, Cosh m(L_-x) 

Jl. f 9 

Lc 0 



dx 



( 17 ) 



9 




mLc 



Tanh 



mLc 



( 18 ) 



9 = 



W 



- T 



SAT 



( 19 ) 



where : 



(T^-T ^^) = average difference between surface 

wall temperature and fluid 
saturation temperature (K) 



Since the unenhanced ends are very long and the thickness 
is very small, only very little heat transfer would take 
place at these ends. Thus, these ends undergo free convection. 



36 



for which the Churchill and Chu correlation as stated in 
[Ref. 10] for the average Nusselt number was assumed: 



Nu, 



1 



0.3 87 Ra, 



1/6 



0.60 + 



I 



I [1 + (0.559/Pr) j 



( 20 ) 



-5 — 12 

10 ^ < Ra^^ < 10 

o 



The average Nusselt number is ; 



Nu, 



h D 



( 21 ) 



where : 



k = thermal conductivity of Freon-113 (W/m.K) 
Solving for h : 



k I 

h = ; 0.60+- 

O 



0.387 

^o 



[1 + (0. 559/Pr) ^/27 j 



( 22 ) 



Pr = - 

a 



(23: 



Ra, 



a B (T - T ) D 
^ ^ W SAT^ o 

V a 



(24) 



B = - — ^ 

p AT 



( 25 ) 
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where : 



Pr 

V 

a 




g 

6 

P 



Prandtl number 

2 

kinematic viscosity (m / s) 

2 

thermal diffusivity (m /s) 
average Rayleigh number 

2 

gravitational acceleration (m/s ) 

volumetric thermal expansion coefficient 

(1/K) 

density of freon-113 (kg/m^) 



Now, substitution of Equations 7, 18 and 24 into Equation 22 
results in: 



h = 



0.387 



hP 1/2 „ 1/6 



g6DlegTanh|j^) ' l/ 



vaL_ ( 



HP 1 1/2 



C 'kA, 



J 



0.60 + 



(26) 



[1 + (0 . 559/Pr) 



Equation 26 is solved for h by iterative technique 
within a range of precision of 0.001. Knowing the value of 
natural-convective heat- transfer coefficient at unenhanced 
ends, the total heat-loss rate is calculated from Equation 
11, and the heat-transfer rate through the enhanced Gewa-T 
surface is obtained by subtracting the total heat-loss rate 
from the total heat-transfer rate: 



^LOSS 






(27) 
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Q 



(28) 



= Qh - 2Q^ 

Finally, the heat flux from the Gewa-T surface at the 
base of the T-shaped fins is: 

q = Q/Ag (29) 

Ag = 7T D 2 L (30) 
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IV. RESULTS AND DISCUSSION 



A. DESIGNATION OF DATA FILES 

A total of seven Gewa-T tubes were tested during this 
investigation. As can be seen from Table 1, the primary 
geometric variables were: fin-tip gap, fin density and 

outside diameter of tubes including fins. 

The general designation for data files is: 

GTNPnn 



where : 

GT indicates Gewa-T surface 

N specifies the number of the tube tested 
(Table 1) 

P indicates the position of the test section 

(i.e., H for horizontal arid V for vertical), 
and 

nn indicates the corresponding cumulative run 
number . 

For example, GT5H28 specifies tube number five tested in the 
horizontal position during run number 28. From Table 1, 
the corresponding specifications for this test section are: 
25.3 mm outside diameter, 23.1 mm diameter at the base of the 
fins, 1.02 fins/mm and inter-fin gap of 0.15 mm. While the 
file name just presented contain data in the very raw form, 
the processed data (i.e., q versus AT) were stored in a 
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separate file. This plot data file could easily be recognized 
by the "P" added at the beginning of the raw data file name. 

In this manner, the plot data file corresponding to the above 
raw data file has the name PGT5H28. 

A total of thirty different data runs were carried out 
during this investigation. These runs break down in the 
following manner: 



Run Number 


Tube Number 


Test Condition 


1-3 


1 


Plain 


- Vertical position 


4-10 


1 


Plain 


- Horizontal 


11 - 21 


1 


With ; 


shrouds and/or Wires 


22 - 23 


2 


Plain 


- Horizontal 


24 - 25 


3 


Plain 


- Horizontal 


26 


4 


Plain 


- Horizontal 


27 - 28 


5 


Plain 


- Horizontal 


29 


6 


Plain 


- Horizontal 


30 


7 


Plain 


- Horizontal 



B. CIRCUMFERENTIAL UNIFORMITY OF HEAT FLUX AROUND HEATER 
ELEMENT 

During the early stages of this investigation, it was 
thought necessary to test if the heater element would provide 
a non-uniform heat flux along the circumferential direction. 
Such a non-uniform heat-flux variation was assumed possible 
based on the heater-element construction. Therefore, tube 1 
was immersed vertically in the R-113 liquid pool and three 
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runs were made. All of these runs showed only a negligible 
variation of the wall thermocouple temperatures; maximum 
variation was ±0.25 K. These runs clearly showed that the 
heater element indeed provided a uniform heat flux along the 
circumferential direction. No detailed data are presented 
in this thesis for the tube in the vertical position. All 
other results presented herein are for the tubes in the 
horizontal position. 

C. LONGITUDINAL UNIFORMITY OF HEAT FLUX ALONG THE TEST 

SECTION 

During most of the data runs with the tubes in the hori- 
zontal position, somewhat unexplainable phenomena were 
observed. As the heat flux was gradually increased, bubble 
formation began simultaneously at the two ends, while the 
center portion of the tube underwent only natural convection. 
Approximately 40 percent of maximum heat flux was necessary 
to obtain nucleate boiling from the entire tube. Since all 
the wall thermocouples are located at midway along the tube, 
the average wall temperature values provided in all data 
runs may be in slight error. 

D. PLOT ANALYSIS FOR PLAIN SURFACE 

Figure 7 represents a typical curve for nucleate pool- 
boiling process of Gewa-T surface in Freon-113. The observed 
behavior of this process is analyzed in this section by 
studying the different sections of the boiling curve. 
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q 



/ 




A-B; Natural Convection 
B-C : Mixed Boiling (Activation of 

Nucleation Sites) 

C-D: Active Boiling 

D-E: High Flux Nucleate Boiling 

E-F; Stable Nucleate Boiling 



Figure 7. Typical Nucleate Pool-Boiling Curve 
for a Gewa-T Surface 
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From point A to point B, a continuous increase in heat 
flux, accompanied by an increase in wall superheat (T^ - 
is observed when heat flux is increased. This region corres- 
ponds to natural-convection process, as also verified during 
experimental runs with the absence of bubbles along both 
enhanced surface and unenhanced ends. 

From point B to point C, a drastic reduction in wall 
superheat was observed while the heat flux continuously 
increased. This region was somewhat unstable as more and 
more nucleation sites became active. This region is known 
as the mixed boiling region, where transition from natural- 
convective boiling to nucleate pool boiling takes place. At 
these heat fluxes, visual observations revealed that the 
bubble formation was quite random. While some portions of 
the tube showed bubbles, other portions showed no bubble 
formation. The wall superheat continued to decrease until 
all nucleation sites became active, as indicated by point D. 

Once all available nucleation sites have been activated, 
the wall superheat starts to increase with increasing heat 
flux as shown by region D to E. When boiling in this region, 
the entire tube showed very high density of bubble formation. 

At this point, it is important to point out that during all 
experimental runs, for the entire range of power settings (7.5 to 
485 W) , little or no nucleate pool boiling on the unenhanced 
surface of the ends was observed (see Figures 8 and 9) for 
the different test sections tested. This observation 
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Figure 8. Photograph of Gewa-T Surface at High Heat Flux 
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Figure 9. Photograph of Gewa-T Surface at Low Heat Flux 




confirmed the validity of the mathematical model assumed for 
data reduction. However, a small rate of bubble formation 
was evident at the tip of the unenhanced fin end; these 
bubbles were in fact generated at the base of the fin (i.e., 
between the Teflon plug and the inner surface of the fin) . 

The vapor bubbles generated in the inter-fin cavities 
are carried up into the channels around the tube due to 
buoyancy forces. As these bubbles travel up, they activate 
more and more nucleation sites as well as they help for the 
removal of larger bubbles from the channels. This phenomenon 
is generally referred to as the "Chimney Effect." 

If the heat flux is too large (see Figure 7) , the rate of 
vapor bubble generation may be too large that the interfin 
channels could be full of a vapor blanket that may disrupt 
the flow of liquid to the hot metal surface where active 
boiling takes places. At this point, further increase in 
heat flux may result in excessive wall temperatures. Under 
these conditions, vapor bubbles were seen ejecting radially 
outward from the test surface. Line D-E represents the boil- 
ing regime just described. It is clear that the slope (i.e., 
heat- transfer coefficient) at E is considerably smaller than 
that at point D. This decrease in heat- transfer coefficient 
is assvimed to be due to the existence of an insulating vapor 
blanket. 

When heat flux is gradually decreased, the curve follows 
a different path (D-F) as shown in Figure 7. This is due 
to the existence of stable nucleation sites that will remain 
active for a wide range of heat fluxes. 
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E. REPEATABILITY WITH PREVIOUS DATA 

Following modifications made, especially to the data 
gathering means (i.e., automatic data-acquisition system and 
newly-written computer programs) , it was necessary to take 
data for tube 1 (this is the same tube tested by Hernandez) 
to check the present data with Hernandez . data under similar 
conditions. Even though a total of 5 different runs were 
made during this stage through debugging of programs and 
system familiarization, only run number 5 will be presented 
in this thesis with a comparison to Hernandez data. 

Figure 10 shows the comparison between the present data 
and the data of Hernandez. A very good agreement is evident 
as the two sets disagree by less than 2 percent for the 
nucleate boiling regime. However, considerable disagreement 
exists in the region of increasing heat flux before total 
nucleate boiling. This is believed to be due to the rather 
random behavior of the boiling process during this regime 
as it will be presented in the next section, even two consecu- 
tive runs for a given tube would not give the same results 
during the unstable boiling regim.e. 

F. REPEATABILITY STUDY OF PRESENT WORK 

Figure 11 shows data from two different runs completed 
on two different days. As discussed earlier, a considerable 
disagreement can be seen before the tube undergoes complete 
nucleate-boiling process. However, for the complete 
nucleate-boiling process, the data from different data runs 
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Figure 10. Reproducibility Test: Hernandez Run No. 3 Versus 

Present Work (GT1H05) 
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Figure 11 . Reproducibility Test: Plain Gewa-T Tube 



agree to within 0.2 percent at a heat flux of 30,000 to 
2 

100,000 W/m . These data runs show considerably good repeata 
bility of this experiment. (Similar agreements were also 
found for all other tubes.) 

Also, these data show a considerable enhancement of 
Gewa-T tube performance over the "Smooth" tube data of 
Hernandez. For example, at AT = 3.1 K, the Gewa-T tube 
gives a 4.65 times larger heat flux than the smooth tube. 

As also discussed by Lepere [Ref. 9] and Hernandez [Ref. 
8] this enhancement in boiling performance is mainly attri- 
buted to the existence of re-entrant cavities in the Gewa-T 
tube . 

G. BOILING PERFORMANCE OF SMOOTH TUBE 

To compare the performance of Gewa-T tubes to smooth- tube 
performance, data of Hernandez were chosen. It is important 
to point out that his smooth tube did not consist of a commer 
cially- available surface. Instead, Hernandez obtained his 
"Smooth" surface by machining fins off of a Gewa-T tube. 

He did not remove fins completely; the remaining portions on 
the t\ibe measured about 0.1 mm. 

Figure 12 shows Hernandez 's data (heat flux as a function 
of wall superheat) for only decreasing heat fluxes. A 
least-squares- fit curve generated for these data points will 
be used in future figures representing data on tube 1. Since 
tubes 2 through 7 have a larger diameter, it was necessary 
to revise Hernandez 's data to include the difference in 
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Smooth-Tube Data from Hernandez [Ref. 7] 
Least-Squares Line for Hernandez Data 
Revised for OD = 23.1 mm Using Equations on 
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Figure 12. Extrapolation to Smooth Tube (23.1 mm OD) from Smooth-Tube 
Data of Hernandez (19.61 mm OD) 



diameter. For this purpose, a semi-empirical correlation 
developed by Cornwell et al. [Ref. 11] is used: 
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Using this relationship, Hernandez 's data were revised 
and a least-squares-f it curve was generated. This curve is 
also plotted in Figure 12. 

H. BOILING PERFORMANCE OF GEWA-T TUBE WITH SHROUDS 

The same four shrouds tested by Hernandez were used 
during this investigation. These shrouds had the following 
geometry : 
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Shroud Number 



Upper Opening 



Lower Opening 



1 


O 

O 


60 


2 


60° 


30 


3 


30° 


30 


4 


O 

o 


8 



All of these shrouds had been designed for tube 1; thus, 
no shrouds were used for tubes 2 through 7 . 

Figures 13, 14, 15, and 16 show data on tube 1 with shroud 
numbers 1, 2, 3 and 4, respectively. As also reported by 
Hernandez [Ref. 8] , all the shrouds improved boiling per- 
formance at low heat fluxes compared to plain Gewa-T tube. 
Table 2 lists performance calculations for tube 1 with various 
shrouds. It can be seen that the 60° x8.5° shroud (number 4) 

produced the maximum improvement of 253 percent in the boiling 

2 

coefficient over smooth tube at a heat flux of 10,000 W/m . 

At this heat flux, the 30° x30° shroud resulted in a 165 

percent improvement in the boiling heat-transfer coefficient. 

2 

At a heat flux of 50,000 W/m , all shrouds resulted in about 
190 percent improvement, except the 30° x30° shroud which 
resulted in a 180 percent improvement. 

2 

At very high heat fluxes, for example at q = 100,000 W/m , 
the 60° x60° shroud showed a maximum performance increase of 
106 percent; the 60° x30° shroud resulted in a 49 percent 
improvement; the 60° x8.5° shroud resulted in a 18 percent 
increase; and the 30° x30° shroud resulted in an 13 percent 
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Figure 13. Effect of 60° x60° Shroud on Boiling Performance 
of Plain Gewa-T Tube 1 
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Figure 14. Effect of 60° x 30° Shroud on Boiling Performance 
of Plain Gewa-T Tube 1 
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Figure 15. Effect of 30° x 30° Shroud on Boiling Performance 
of Plain Gewa-T Tube 1 
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Figure 16. Effect of 60° x8.5° Shroud on Plain Gewa-T Tube 
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improvement. It can be seen that the values just presented 
are somewhat proportional to the summation of the two shroud 
angles. This behavior indicates that at very high heat 
fluxes, too small a shroud opening can obstruct the outflow 
of vapor and/or the inflow of liquid. 

However, it must be noted that at all practical heat 
fluxes (10,000 to 80,000 W/m^) , all the shrouds resulted in 
a considerable enhancement in the boiling performance. The 
60° x 8.5° shroud appears to be the best selection for all 
practical heat fluxes. 

I. PERFORMANCE OF GEWA-T TUBE WITH WIRE WRAPS IN INTER-FIN GAP 

A series of runs were made with two to five 0.1 ram-diameter, 
copper wire wraps in the inter-fin gap of tube 1. Figures 
17, 18, 19 and 20 show data with wire wraps of 2, 3, 4 and 
5, respectively. It can be seen that all of these curves 
show increases in boiling performance at all wall superheats 
for decreasing heat flux values. This enhancement is most 
probably due to the increase in the nucleation sites within 
the inter-fin cavities. However, there is a small reduction 
in performance in the natural convection process probably 
owing to that wires could act as an insulating material due 
to the thermal contact resistance between the wires and the 
inter-fin area. 

Table 2 also shows that the use of three wires resulted 
in the best performance increase. This combination resulted 
in a 341 percent improvement in the boiling heat-transf er 
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Figure 17. Effect of Two Wires Wrapped in the Inter-fin 
Cavity of Tube 1 on Boiling Performance 
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Figure 18. Effect of Three Wires Wrapped in the Inter-fin 
Cavity of Tube 1 on Boiling Performance 
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Figure 19. Effect of Four Wires Wrapped in the Inter-fin 
Cavity of Tube 1 on Boiling Performance 







64 



Figure 20. Effect of Five Wires Wrapped in the Inter-fin 
Cavity of Tube 1 on Boiling Performance 




coefficient over the smooth tube value at a wall superheat 
of 1.4 K. 

J. BOILING PERFORMANCE OF TUBE 1 WITH 60° x8.5° SHROUD AND 

THREE WIRES 

Data were also taken to combine the best shroud and the 
best number of wire wraps in the inter-fin channel. For 
this purpose, three wire wraps were added to the inter-fin 
channel and data were taken with the 60° x 8.5° shroud. As 
Figure 21 and Table 2 show, this combination resulted in a 
performance yet superior to the use of three wires or the use 
of 60° x8.5° shroud separately. Table 2 shows that this 
arrangement resulted in 481 percent increase in the boiling 
heat-transfer coefficient over the smooth-tube value at a 
heat flux of 10,000 W/m^ . 

K. EFFECTS OF FIN-TIP GAP AND FIN DENSITY OF GEWA-T 
. TUBES ON BOILING PERFORMANCE 

As described earlier, six Gewa-T surfaces were tested 
during this stage of the investigation; tubes 2, 3, and 4 
had a fin density of 0.75 fins/mm, and tubes 5, 6 and 7 
with a fin density of 1.02 fins/mm, and each set of tubes 
with fin-tip gaps of 0.15, 0.25 and 0.33 (see Table 1). 

Micrographs were obtained for all enhanced surfaces on 
a Scanning Electron Microscope with a magnification of 50X. 
Figure 5 shows micrographs representing tubes 2 and 5. 
Computed surface areas and area ratios (with respect to 
smooth tube) are listed in Table 3. 
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Figure 21. Effect of Three Wires Wrapped in the Inter-fin Cavity 
with 60° x 8.5° Shroud on Gewa-T Tube 1 on Boiling 



TABLE 3 



Computed Surface Areas and Area Ratios 



Tube 

1 

2 

3 

4 

5 

6 
7 



Cavity 

Surface 


and Fin 2 
Area (m ) 


Total Surface Area 
(m^) 


^F'^'^SMOOTH 


240.47 


X lo"^ 


8.818 


X lO"^ 


— 


288.01 


X lO"^ 


10.944 


X lo"^ 


2.74 


305.34 


X lO"^ 


11.603 


X lO"^ 


2.91 


285.07 


X lO"^ 


10.833 


X lO"^ 


2 . 72 


256.04 


X 10“^ 


13 . 314 


X lO"^ 


3 . 60 


248.61 


xio'^ 


12.928 


X lO"^ 


3.50 


267.91 


X lo"^ 


13.931 


X lO"^ 


3 .77 



It can be seen that tubes with a higher fin density 
(1.02 fins/mm) result in larger area ratios than the tubes 
with lower fin density (0.75 fins/mm) . However, for a given 
fin density, the variation of the area ratio with fin-tip 
gap is not conclusive. 

Figures 22 through 27 show data taken on tubes 2 through 
7, respectively. These results are also listed in Table 2. 
This table shows that tube 6 (with a fin density of 1.02 
fins/mm and a fin-tip gap of 0.25 mm) resulted in the highest 
heat flux at all values of wall superheat. 

To investigate the effects of fin density and fin-tip 
gap on boiling heat-transfer performance, the results pre- 
sented in Figures 22 through 27 (which showed variations of 
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Figure 22 . Boiling Performance of Gewa-T Tube 
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Figure 23. Boiling Performance of Gewa-T Tube 
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Figure 24. Boiling Performance of Gewa-T Tube 
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Figure 25. Boiling Performance of Gewa-T Tube 
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Figure 26. Boiling Performance of Gewa-T Tube 
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Figure 27. Boiling Performance of Gewa-T Tube 



heat flux as a function of wall superheat) were re-plotted 
on a different basis as shown in Figure 28. This figure shows 
the variations of the enhancement of heat-transf er coefficient, 
over the smooth-tube value, with wall superheat for all six 
tubes. It is evident that tube 6, which has a fin density 
of 1.02 fins/mm and a fin-tip gap of 0.25 mm results in the 
maximum enhancement at all wall superheat values. 

Tube 5 shows poorer performance than tube 6. This poorer 
performance is believed to be due to the fact that the smaller 
fin-tip gap disrupts the flow of liquid into the channels 
and/or the vapor flow out of the channels. Tube 7 (which has 
the largest fin-tip gap of 0.35 mm) shows the poorest per- 
formance among the set of these tubes with a fin density of 
1.02 fins/mm. As visually observed during experimental 
runs, this tube ejected a considerable portion of vapor 
bubbles radially along tube periphery. This premature 
departure of vapor bubbles results in a lower performance 
by producing a thiclcer liquid film at the hot surface. 

In comparison with the results for tubes 5, 6 and 7, the 
curves representing tubes 2, 3 and 4 show that for a given 
fin-tip gap, the performance has decreased with decrease in 
fin density. This observation is believed to be due to the 
presence of a greater number of channels (for a given length 
of tube) in the higher-fin-density tubes. Again, for a fin 
density of 0.75 fins/mm, the 0.25 mm fin-tip gap (tube 3) 
resulted in the best boiling performance. The same 
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Figure 28. Enhancement Ratios of Gewa-T Tubes over Smooth 
Tubes Based on Nominal Area 



explanations given above should be valid for the poorer per- 
formance of tubes 2 and 4. 

It can be seen that all tubes show smaller enhancement 
ratios at low wall superheats, and the above-mentioned 
differences in performance are more pronounced. On the 
other hand, all the tubes show greater enhancements at larger 
wall superheats and they seem to approach a common enhance- 
ment ratio. The increase in enhancement is believed to be 
due to the increased chimney effect, which helps to remove 
vapor bubbles from channels resulting in a thinner liquid 
film at the hot boiling surface. When heat flux is further 
increased, vapor can provide a blanket around the hot 
surface, which may eventually result in a critical heat 
flux. It appears that the critical heat flux is not strongly 
dependent on the fin density or fin-tip gap. 

The results shown in Figure 28 were also replotted as 
shown in Figure 29 to further study the performance of these 
six tubes. This figure shows the variation in enhancement 
ratio obtained beyond the area increase (i.e., the ratio of 
heat-transfer enhancement/area ratio) as a function of wall 
superheat. It is evident that obviously all tubes resulted 
in enhancement ratios beyond the area increase. 

In summary, a fin density of 1.02 fins/mm resulted in 
performance superior to a fin density of 0.75 fins/mm. Also, 
in both cases, a fin-tip gap of 0.25 mm resulted in the best 
performance . 
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Figure 29. Enhancement Ratios of Gewa-T Tubes over 
Smooth Tube Based on Actual Area 



V. CONCLUSIONS 



(1) Experimental observations revealed a non-uniform 
temperature distribution along test section. The bubble 
generation was seen to start at the two ends and migrate into 
the center as the heat flux was increased. No clear explana- 
tion is possible for this temperature distribution. 

(2) Results of the present investigation showed very 
good repeatability with the results of Hernandez obtained 
under similar conditions. 

(3) The use of shrouds improved the boiling heat-transf er 

coefficient of Gewa-T surface at all practical heat fluxes. 

However, a considerable reduction in heat-transfer coefficient 

2 

was seen at very high heat fluxes (> 100,000 W/m ) . Among 
the four shrouds tested, the shroud with 60° x8.5° openings 
resulted in the best performance verifying Hernandez's 
results. This shroud resulted in a 253 percent increase in 
boiling heat-transf er coefficient, over smooth tube, at a 
wall superheat of 1.8 K. 

(4) The use of 0.1 mm-diameter wires (2 to 5) improved 
boiling heat-transfer coefficient at all practical heat 
fluxes. The use of 3 wires resulted in the best performance 
with a 341 percent increase in the heat-transf er coefficient 
(over the smooth-tube value) at a wall superheat of 1.4 K. 

(5) Further improvement in boiling performance was 
obtained by using the 60° x8.5° shroud with three wire wraps. 
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This combination resulted in a 481 percent increase in boil- 
ing heat-transfer coefficient at a wall superheat of 1.1 K. 

(6) The use of wires appears to be superior to the use 
of more expensive shrouds. 

(7) Among the fin-tip gaps of 0.15, 0.25 and 0.35 mm 
tested, the 0.25 mm gap resulted in the best boiling 
performance . 

(8) The Gewa-T tubes with a fin density of 1.02 fins/mm 
performed superior to a fin density of 0.75 fins/mm. At a 
wall superheat of 3.5 K, the lower-fin-density tube increased 
the boiling heat- transfer coefficient by 80 percent, while 
the higher-fin-density tube increased by 103 percent. 
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VI . RECOMMENDATIONS 



(1) Investigate the reason for the existence of a 
longitudinal temperature profile and take appropriate 
steps to alleviate this problem. 

(2) Study the possibility of adding roughness in the 
inter-fin cavity of Gewa-T tubes. This operation should 
be performed before the rolling process to obtain the T- 
shaped fins. 

(3) Repeat runs made during this investigation with 
water as the boiling medium. 
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APPENDIX A 



' 000 ! 
1010 ! 
1 P20’ 
1 030! 
10^0 
1050 
I ObO 
1070 
1080 
1090 
1 100 
1110 
I 120 
1130 
I MO 
1150 
1 160 
1170 
I 1B0 
1 190 
1200 
1210 
1220 
1230 
12^0 
1250 
1260 
1270 
1280 
1290 
1300 
1310 
1320 
1330 
13^0 
1350 
1360 
1370 
13B0 
1390 
MOO 
M10 
M20 
M30 
M^O 
M50 
M60 
M70 
M80 
M90 
1500 
1510 
1520 
1530 
15^0 
1550 
1560 
1570 
1580 
1590 



DATA REDUCTION PROGRAM 



-ILE NAME: DRP2 

DATE: May ^3. 198^ 

REVISED: Juiy 20. !984 



BEEP 

PRINTER IS 1 

^^RINT USING . ‘*’*Seiect option: 



PRINT USING ”6X.*”‘0 Taking data or re-process ing previous data 

PRINT USING '*6X*"*’1 Plotting data 

INPUT Idp 

IF ldp=0 THEN CALL Main 
IF Idp-1 THEN CALL Plot 
END 

SUB Main 
COM /Cc/ C(7) 

DIM Emf <9) ,T(9) 

DATA 0.10086091 . 25727 . 94369 , -7873^5 . *3295 . 78025595 . 31 
DATA -9247486589, 6. 97688E+1 1 , -2 .661 92E+1 3 . 3 . 94Q78E^1 4 
READ C(-) 

BEEP 

INPUT '’SELECT TUBE OD (0 = 0.75 in, 1 = 1.0 in)'Mtd 
IF ltd=0 THEN 

DM.0159 ! Diameter at thermocouple positions 

D2=.0193 ? Diameter of test section to the base of fins 

Di=.0183 ! Inside diameter of unenhanced ends 



Do=.0190 



Outside diameter of unenhanced ends 



L=.0490 

Lu=.0319 

ELSE 

DM.01786 
D2=.0231 
Di=. 02032 
Do=. 02222 



Length of enhanced surface 

Length of unenhanced surface at the ends 



L=.0508 
Lu=. 03175 
END IF 

Lc=Lu-*-( Do-D L ) /2 



A=PM(Do 2-Dr 2)/4 
P=PMDo 



Kcu=385 ! Thermal conductivity of Copper 

PRINTER IS 701 
CLEAR 709 
BEEP 

INPUT '‘ENTER MONTH. DATE AND TIME ( MM : DD : HH : MM : SS ) ". DateS 
OUTPUT 709 :"TD": Dates 
OUTPUT 709:"TD“ 

ENTER 709:DateS 

PRINT '* Month, date and time :":DateS 

PRINT 

PRINT USING '*1 OX . "’’NOTE : Program name : BOIL 

BEEP 

INPUT "ENTER DISK NUMBER". Dn 

PRINT USING "16X.""Disk number = ”".ZZ":Dn 

BEEP 



INPUT "ENTER INPUT MODE ( 0 = 3054A , ! =FILE ) '’ , Im 
Im= Im+ 1 
IF Im=l THEN 
BEEP 

INPUT "GIVE A NAME FOR THE RAW DATA FILE" . D._^ i leS 
PRINT USING "16X.""Ne«ii +"iie name: , 1 4A" : D_f l ieS 
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ISOO 

ISIO 

1620 

1630 

1640 

1650 

1660 

1670 

1680 

1690 

1700 

1710 

1720 

1 730 

1740 

1750 

1 760 

1770 

1780 

1790 

1800 

1810 

1820 

1830 

1840 

1850 

I860 

1870 

1880 

1890 

1900 

1910 

1920 

1930 

1940 

1950 

I960 

1970 

1980 

1990 

2000 

2010 

2020 

2030 

2040 

2050 

2060 

2070 

2080 

2090 

2100 

2110 

2120 

2130 

2140 

2150 

2^60 

21 70 

2180 

2190 



BEEP 

INPUT ’’ENTER PLUID < Q = f^PEOM- I 1 3 . 1 ='^ATER ) ” . Cf 
Cf-Cf+1 

CREATE BOAT D_riUS.2G 
ASSIGN '?File TO D.fileS 
OUTPUT «Fiie:DateS.Cf 
ELSE 
BEEP 

INPUT ”GIUE THE NAME OF THE EXISTING DATA FILE” . D_f i ieS 

BEEP 

lhc = 0 

INPUT ’’ENTER 1 IF USING FILE NUMBER < 25’Mhc 
PRINT USING ”16X.””Qld file najne : ”’M 4A” ; i leS 
BEEP 

INPUT "ENTER THE NUMBER OF RUNS STORED”. Nran 
ASSIGN '^File TO D_fileS 
ENTER ?File:DoldS.Cf 

PRINT USING "16X.””This oatpat for file created at . 1 4A” : Do ld$ 

END IF 
BEEP 

INPUT "GIVE A NAME FQR PLOT FILE” . P_f i 1 eS 
CREATE BDAT P_fileS.5 
ASSIGN apiot TO P_fileS 

IF Cf=l then print USING ” 1 6X . ””F Imd : Freon-113 

IF Cf=2 THEN PRINT USING "16X."”Flaid : Water 

BEEP 

INPUT "ENTER OUTPUT VERSION < 1 =SHORT . 0=LONG ) ” . lov 

IF Iov-0 THEN Iov=2 

BEEP 

INPUT "ENTER NUMBER OF DEFFECTIVE TCS < 0 = DE.FAULT ) ” . Idtc 

IF ldtc>0 THEN 

BEEP 

IF Idtc=1 THEM 

INPUT "ENTER DEPi^ECTIVE TC LOCATION” . Ldtc 1 

PRINT USING ”16X.””TC is deffective at location ””,D”:Ldtcl 
END IF 

IF Idtc-2 THEN 

INPUT "ENTER DEFFECTIVE TC LOCATIONS” . Ldtc I . Ld tc2 

PRINT USING "16X."”TC are deffective at locations . D , 4X . D” :Ldtc1 . Ldtc2 

END IF 

END IF 

IF Im=I THEN 

BEEP 

INPUT "ENTER BAROMETRIC PRESSURE < inHg) " . Patm 
Patm=Patm+2* 1500/13500 
OUTPUT ^F lie: Patm 
ELSE 

ENTER ^File:Patm 

IF Ihc = 1 THEN Patm=Patm-*-2*-1 500/1 3500 
END IF 

Patm* Patm *25 .4*14. 696/760 
Tsat=FNTsat<Cf .Patrn) 

PRINT USING " 1 6X . ""Sat lirat ion temperature = ””.3D.DD."” (Deg C) : Tsat 

J=1 
Sx*0 
Sy-0 
S X s = 0 
S X V = 0 
BEEP 

INPUT "HIT CONTINUE WHEN READY”. Ok 
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2200 Repeat: I 

2210 IF Lti=! them 

2220 OUTPUT 709:”PiR PF29 AL38 WS" 

2230 FOR 1=0 TO 9 
22^0 OUTPUT 709;”AS SA“ 

2250 Sam=0 
2260 FOR Ji=1 TO 20 
2270 ENTER 709;E 
2280 SLLm = S'Jjn*^E 
2290 NEXT Ji 
2300 Emf (I)=Slud/20 
2310 NEXT I 
2320 BEEP 

2330 INPUT -ENTER HEATER VOLTAGE. Vh (V)’\Vh 
2340 BEEP 

2350 INPUT -ENTER RESISTOR VOLTAGE. Vr (V)-.Vr 
2360 ELSE 

2370 ENTER ^F i le : Emf ( - ) , Vr . Vh 

2380 END IF 

2390 Tua=0 

2400 FOR 1=0 TO 9 

2410 IF ldtc>0 THEN 

2420 IF (I + 1)=Ldtc1 OR (I + D=Ldtc2 THEN 

2430 T(I)=99.99 

2440 GOTO 2490 

2450 END IF 

2460 END IF 

24 70 T(I)=FNTvysv(Emf (I) ) 

2480 IF I<8 THEN Tua=T.ja+T( I) 

2490 NEXT I 

2500 Tua=Tiiia/ < 8- Id tc ) 

2510 Tld=T(9) 

2520 Tv=T(8) 

2530 Rho=FNRho < Cf . Tsat ) 

2540 Mti=FNMu(Cf .Tsat) 

2550 K=FMK(Cf .Tsat) 

2560 Cp=FNCp < Cf , Tsat ) 

2570 Beta=FNBeta( Cf . Tsat ) 

2580 Ni=MtjL/Rho 
2590 A lpha=K/ < Rho*Cp ) 

2600 Pr=Ni/Alpha 
2610 Amp=Vr/2.031 

2620 0=Vh*Amp 

2630 T.j=Tua-Q*L0G(D2/Dl ) /(2 *PI-Kcul*L) 

2640 The tab= Till- Tsat 
2650 IF ThetaD<0 THEN 
2660 BEEP 

2670 OISP -TWALL<TSAT (HIT CONTINUE WHEN READY)- 

2680 PAUSE 

2690 GOTO 2210 

2700 END IF 

2710 Hbar=l90 

2720 Fe= ( Hbar*<^/ ( K cll*A ) ) .5*Lc 
2730 Tanh=FNTanh(Fe) 

2740 Theta=ThetaD*Tann/Fe 

2750 Xx= ( 9 . 31 ♦Beta^The tab*Do 3*Tanh/ ( Fe^Mi -A Ipha) ) .166667 
2760 Yy = (U( .559/Pr) (9/16)> (8/27) 

2770 Hbarc = K/Do*( .6+.387-Xx/Yy ) 2 
2780 IF ABS ( ( Hbar-Hbarc) /Hbarc ) > . 00 1 THEM 
2790 Hbar= ( Hbar+Hbarc) ^ . 5 
2800 GOTO 2720 
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2310 END IF 

2320 Q i = ( ) . S^Thetab^Tann 

2830 Qc=Q-2-Qi 
2340 As=PI-D 2*L 
2850 Qdp=Qc/As 

2360 IF Iov=i THEN PRINT USING •’ ! OX . DD , 2X . 3D . DD . 2X . 2( Z , 4DE . 2X ) . DD , DD*’ : J . Tua. Hba 
r . Odp . The tab 
,2370 II- Iov = ^ THEN 

2880 PRINT USING “ ! OX . •'-Data set namber = ‘”’.DD”:.j 

2890 PRINT USING " 1 OX . ‘"’Pos 1 1 ion nuunber : 1 2 3 4 5 6 

1 8 

2900 PRINT USING ” 1 OX . •’•'Temperature : ’'•' ,3( DD .DD , 1 X) " : T( 0 ) , TO ) . T< 2) . T( 3) . T 
(4) . T(5) . T(6) . T( 7) 

2910 PRINT USING " 1 OX , ‘"*T1 igd Tvapr Tua Hbar Qdp Thetab’"'” 

2920 PRINT USING ”9X ,3(3D . DD .2X) .2(Z.4DE.2X) . 2D . DD” : T Id . Tv . Twa. Hbar . OdP . Thetab 

2930 END IF 

2940 IF Im-1 THEN 

2950 BEEP 

2960 INPUT ’'OK TO STORE THIS DATA SET < 1 - Y , 0=N ) . Ok 

2970 END IF 

2980 IF Gk=l OR Im-2 THEN J=J+1 

2990 IF Ok = l AND Im-1 THEN OUTPUT '^F i ie : Emf ( * ) . Vr . Uh 
3000 IF Im=2 OR Ok = 1 THEN OUTPUT «P lo t : Odp , The tab 
3010 IF Im=1 THEN 
3020 BEEP 

3030 INPUT ”NILL THERE BE ANOTHER RUN ( t =Y . O’N ) . Go_on 

3040 Nrun=J 

3050 IF Go_on01 THEN 3110 
3060 WAIT 300 

3070 IF Go_on= I THEN Repeat 
3080 ELSE 

3090 IF J<Nrun+1 i HEN Repeat 

3100 END IF 

3110 IF Im=i THEN 

3120 BEEP 

3130 PRINT 

3140 PRINT USING ”I0X2'”N0TE: ’"'.ZZ.’”' data runs -jere stored in file ’"'.10A”:J- 

1 ,D^f lies 

3’50 END IF 

3160 BEEP 

3(70 PRINT 

3(30 PRINT USING ” 1 OX . ’’"NOTE : “”.ZZ.‘”' X-Y oairs uere stored in plot data file 
. lOA” : J-1 ,P_fileS 
3190 ASSIGN iiFile TO - 
3200 ASSIGN t^Plot TO - 
3210 BEEP 

3220 INPUT ’'LIKE TO PLOT DATA ( 1 = Y , 0 = N ) ?" . Ok 

3230 IF 0k=l THEN 

3240 CALL Plot 

3250 END IF 

3260 SUBEND 

3270 DEF FMMuCCf .T) 

3280 IF Cf=l THEN 

3290! RANGE: 28 TO 50 DEG C 

3300 Mu=8. 96298 19E-4-T-( 1 , 1 094609E-5-T-5 . 566829E-8 ) 

3310 END IF 

3320 IF Cf=2 THEN 

3330 A=247.3/(T+1 33,15) 

3340 Nu=2.4E-5-l0 A 

3350 END IF 
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336G RETURN Nu 

3370 FNEND 

3380 DEF rNCp(Cf.n 

3390 IF Cf=1 THEN 

3400! RANGE: 0 TO 30 DEG C 

3410 Cp=9.2507273E-l+T-<9.3400433E-4+l . 7207792E-6-T) 

3420 END IF 
3430 IF Cf=2 THEN 

3^140 Cp=4.21120858-T-(2.26826E*3-T-(4.^2361E-5+2.7t428E-7-T) ) 

3450 END IF 

3460 RETURN Cp*1000 

3470 FNEND 

3480 DEF FNRho(Cf.T) 

3490 IF Cf=1 THEN 

3500! RANGE: 30 TO 80 DEG C 

3510 Ro=1 . 6207479E+3-T*< 2. 2186346+1*2. 3578291 E- 3) 

3520 END IF 
3530 IF Cf=2 THEN 

3540 Ro=999.52946+T-< . 0 1 269-T-( 5 . 4825 1 3E-3-T- 1 .234t47E-5>) 

3550 END IF 

3560 RETURN Ro 

3570 FNEND 

3580 DEF FNPr(Cf.T) 

3590 Pr-FNCp(Cf .T)*FNMu<Cf .T)/FNK(Cf .T> 

3600 RETURN Pr 

3610 FNEND 

3620 DEF FNK(Cf,T) 

3630 IF Cf=l THEN 

3640! RANGE: 30 TO 90 DEG C 

3650 K=8.2095238E-2-T*(2,2214286E-4+T-2.3809524E-8) 

3660 END IF 

3670 IF Cf=2 THEN 

3680 X-<T+273. i5>/273. 15 

3690 K = -.92247+X^(2.8395-X»( I .8007-X*( . 52577- . 07344*X ) > > 

3700 * END IF 
3/10 RETURN K 
3720 FNEND 
3730 DEF FNTanh(X) 

3740 P=EXP<X) 

3750 0=EXP(-X) 

3760 Tanh=(P-Q>/(P+rg) 

3770 RETURN Tanh 

3780 FNEND 

3790 DEF FNTvsv(U) 

3800 COM /Cc/ C(7) 

3810 T=C(0) 

3820 FOR 1-1 TO 7 
3830 T-T + Cd)-^; I 
3840 NEXT I 

3850 T-T+8 .34643E- I -T*( 2 . 4547896E-2-T* 1 . 3722583E-^> 

3860 RETURN ^ 

3870 ^NEND 

3830 DEF FNBeta(Cf.T) 

3890 Rop=FNRho(Cf ,T+. 1 ) 

3900 Rom=FNRho(Cf d-. 1 > 

3910 Be ta=-2/ < Rop+Rom ) ♦ < Rop-Rom) / .2 
3920 RETURN Beta 
3930 FNEND 

3940 DEF FNTsat(Cf.P) 

3950 Tl-10 

2960 Th=110 
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3370 ia=(Tl^Th>-.5 

3380 Pc=FNPsar(Cf .Ta) 

3990 IF ABS( (P'Pc)/P)> .0001 THEN 

^000 IF Pc>P THEN Th=Ta 

4010 IF Pc<P ^HEN Tl=Ta 

4020 GOTO 3970 

4030 END IF 

4040 RETURN Ta 

4050 FNEND 

4060 DEF FNPsat<Cf .Tc) 

4061 IF Cf=1 THEN 
4070 T=Tc-l .8+32+459.6 

4080 P=10 M33. 0655-4330. 38/T-9.2635*LGT( T) +2. 0539E-3*T> 

4081 END IF 

4082 IF Cf=2 THEN 

4083 P=Pvstu(Tc) 

4085 END IF 
4090 RETURN P 
4100 FNEND 

4110 SUB Plot 
4120 DIM C(9) 

4130 INTEGER Ii 
4140 PRINTER IS 1 
4150 Idv=1 
4(60 BEEP 

4170 INPUT '‘LIKE DEFAULT UALUES FOR PLOT ( 1 = Y . 0 = N > ^’Mdv 
4180 BEEP 

4(90 PRINT USING "4X , '"‘Select Option:'""* 

4200 PRINT USING "4X.""0 q versus delta-T 

4210 PRINT USING '*4X,*"*1 h versus delca-T 

4220 PRINT USING "4X.’"'2 h versus q 

4230 PRINT USING "4X,'"*3 h-ratio versus delta-T 

4240 INPUT Opo 

4250 IF 0po=3 THEN 
4260 BEEP 

4270 INPUT "SELECT TUBE DIAMETER (0=. 75. 1=1.0 IN)". ltd 

4280 END IF 

4290 PRINTER IS 705 

4300 IF IdvOI THEM 

4310 BEEP 

4320 INPUT "ENTER NUMBER OF CYCLES FOR X-AXIS". Cx 
4330 BEEP 

4340 INPUT "ENTER NUMBER OF CYCLES FOR Y-AXIS". Cy 
4350 BEEP 

4360 INPUT "ENTER MIN X-^;ALUE (MULTIPLE OF IO)",XfTun 
4370 BEEP 

4380 INPUT "ENTER MIN Y-UALUE (MULTIPLE OF IO)".Yfrun 

4390 ELSE 

4400 Cy=3 

4410 IF 0do=0 them 

4420 Cx=2 

4430 XiTun=1 

4440 Ymin=(000 

4450 END IF' 

4460 IF 0po=( then 
4470 Cx=2 
4480 Xmin=i 
4490 Ymin=100 
4500 END IF 
4510 IF 0po=2 THEM 
4520 Cx=3 
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4530 

4540 

4550 

4560 

4570 

4580 

4590 

4600 

4610 

4620 

4630 

4640 

4650 

4660 

4670 

4680 

4690 

4700 

4710 

4720 

4730 

4740 

4750 

4760 

4770 

4780 

4790 

4800 

4310 

4820 

4830 

4840 

4850 

4860 

4870 

4880 

4890 

4900 

4910 

4920 

4930 

4940 

4950 

4960 

4970 

4980 

4990 

5000 

5010 

5020 

5030 

5040 

5050 

5060 

5070 

5080 

5090 

5^00 

5110 

5»20 

5130 



Xmin = '! 000 

Y mi n = ’ 0 0 
END IF 

IF Gpo=3 THEN 
Cx*2 
Cy = 3 
Xm i n= 1 

Y m 1 n = 1 
END IF 
END IF 
BEEP 

PRINT "INiSPliIP 2300,1800,3300,6800;” 

PRINT ”SC 0. 100,0, 100;TL 2,0;" 

Sfx=l00/Cx 

Sfy=100/Cy 

PRINT "PU 0,0 PD" 

Nn=9 

FOR I=I TO Cx^-1 
Xat-Xm in*1 0" < 1“ 1 ) 

IF I=Cx+l THEN Nn=l 
FOR J=1 TO Mn 

IF J=1 THEN PRINT "TL 2 0" 

IF J-2 THEN PRINT "TL 1 0" 

Xa=Xat^J 

X=LGT(Xa/Xmin)*Sf X 
PRINT "PA":X,",0: XT:" 

NEXT J 
NEXT I 

PRINT "PA I00.0:PU:" 

PRINT "PU PA 0.0 PD" 

Nn = 9 

FOR I=I TO Cy+1 
Yat= Ym 1 n*1 0" < I - 1 ) 

IF I=Cy+l THEM Nn=l 
FOR J=1 TO Nn 

IF J=1 then PRINT "TL 2 0" 

IF J=2 THEN PRINT "TL 1 0" 

Ya=Yat-J 

Y=LGT< Ya/Ymin)*Sfv 
PRINT "PA 0,":Y,"YT" 

NEXT J 
NEXT I 

PRINT "PA 0.100 TL 0 2" 

Nn=9 

FOR I=I TO Cx+1 
Xat = Xm in*l 0" ( 1-1 > 

IF I=Cx-l then Nn=l 
FOR J=1 TO Nn 

IF J=1 THEN PRINT "TL 0 2" 

IF J>1 THEN PRINT "TL 0 1" 

Xa=Xat-J 

X=LGT (Xa/Xmin)*STx 
PRINT "PA";X,", 100: XT" 

NEXT J 
NEXT I 

PRINT "PA 100,100 PU PA 100,0 PD" 

Nn = 9 

FOR 1=1 TO Cy+1 
Yat = Yrrun*1 O' ( 1-1 ) 

IF I=Cy-^l THEN Mn-1 
FOR J*1 TO Nn 
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5140 IF J=l "HEN PRINT "TL 0 2** 

5150 IF J>1 THEN PRINT •‘TL 0 I" 

5160 ra=Yat-J 

51 70 Y = LGT( Ya/YfT>in)-Sfy 

5130 PRINT '’PD PA I 00 Y . “YT” 

5190 NEXT J 

5200 NEXT I 

5210 PRINT *'PA 100 . lUO PU” 

5220 PRINT ’’PA 0.-2 SR 1 .5,2” 

5230 I i =LGT < Xm in ) 

5240 FOR 1 = 1 TO Cx-^1 
5250 Xa=Xir>in-IO'(I-1 ) 

5260 X=LGT(Xa/Xmin)-Sf X 
5270 PRINT •*PA”:X,”.0:'‘ 

5280 IF li> = 0 THEN PRINT '’CP -2 . -2 : LB 1 0 : PR -2 ,2 ; LB” ; 1 1 ; 

5290 IF li<0 THEN PRINT ”CP -2 . -2 ; LB1 0 : PR 0 , 2 ; LB” ; I i : 

5300 Ii=Ii+1 

5310 NEXT I 

5320 PRINT ”PU PA 0.0” 

5330 Ii = LGT(YiDin> 

5340 Y10=I0 

5350 FOR 1=1 TO Cy+1 

5360 Ya=Yn)in*10<I-1) 

5370 Y=LGT(Ya/Yi7>in)*$fy 
5380 PRINT ”PA 0,”:Y,”” 

5390 PRINT ”CP -4 . - . 25 : LBl 0 : PR -2 . 2 ; LB” : I i : 

5400 Ii = Ii^-1 
5410 NEXT I 
5420 IF IdvOI THEN 
5430 BEEP 

5440 INPUT ”ENTER X-LABEL” . X Iat>eI3 
5450 BEEP 

5460 INPUT -ENTER Y-LABEL” , Y laoe IS 

5470 ELSE 

5480 IF Opo=0 THEN 

5490 XlabelS = **T.j-Tsat (K)” 

5500 YlabelS=”q (W/m 2)” 

5510 END IF 

5520 IF Gpo=^ THEN 

5530 X 1 abe 1S==” Tu- Tsat (K)” 

5540 Ylaber5=”h (H/m 2.K)” 

5550 END IF 

5560 IF Gpo=2 THEN 

5570 Xlabel$=”a (W/m 2>” 

5580 YlabelS="h (W/m 2.K)” 

5590 END IF 

5600 IF Cpo=3 THEN 

5610 X labe lS = ”Tij-Tsat (K)” 

5620 Y labe lS=”h ( enh ) /h( smooth ) ” 

5630 END IF 
5640 END IF 

5650 PRINT ”SR !.5.2:PU PA 50. -16 CP” : -l.EN< X labe IS) /2 : ”0 : LB” : X labe IS 
5660 PRINT ”PA -14.50 CP 0 .”: -LEN< Y labe IS ) 72-5/6 : ”DI 0 . 1 ; LB” : Y labe IS 
5670 PRINT ”CP 0.0 DI” 

5680 Repeat:! 

5690 BEEP 

5700 INPUT ”WANT TO PLOT DATA FROM A FILE ( ? =Y . 0=M) P” .Ok 
5710 IF 0k=1 THEN 
5720 BEEP 

5730 INPUT “ENTER THE NAME OF 'HE DATA FILE” . i leS 
5740 ASSIGN '?<File TO D^fileS 
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5750 

5760 

5770 

5780 

5790 

5800 

5310 

5820 

5830 

58^0 

5850 

5860 

5870 

5880 

5890 

5900 

5910 

5920 

5930 

59^0 

5950 

5960 

5970 

5980 

5990 

6000 

6010 

6020 

6030 

6040 

6050 

6060 

5070 

6080 

6090 

6100 

6110 

6120 

6130 

6140 

6150 

6160 

6170 

6180 

6190 

6200 

6210 

6220 

6230 

6240 

6250 

6260 

6270 

6280 

6290 

6300 

6310 

6320 

6330 

6340 

6350 



BEEP 

BEEP 

INPUT “’ENTER THE BEGINNING RUN NUflBER’Mid 
BEEP 

INPUT ’’ENTER THE NUMBER OF X-Y PAIRS STORED” . Npai rs 
BEEP 

PRINTER IS 1 



USING 


"4X," 


"Se 


lec t a symDo i .* 


USING 


"4X." 


" 1 


Star 2 PI as Sign 


USING 


"4X." 


"3 


Circle 4 pqaare 


USING 


"4X," 


"5 


Rombus 


USING 


"4X," 


"6 


R ight-s ide“op triangle 


USING 


"4X," 


"7 


Up-s ide-douin triangle 



SM*” 
SM + ” 
SMo” 



INPUT Sym 
PRINTER IS 705 
PRINT ”PU DI” 

IF Sym=1 THEN PRINT 
IF Sym=2 THEM PRINT 
IF Sym=3 THEN PRINT 
IF Md>1 THEN 
FOR 1=1 TO (Md-1 ) 

ENTER «File;Ya,Xa 
NEXT I 
END IF 

FOR 1=1 TO Npairs 

ENTER ^File;Ya,Xa 

IF 0po=1 then Ya=Ya/Xa 

IF 0po=2 THEN 

Q = Ya 

Ya=Ya/Xa 

Xa=0 

END IF 

IF 0po=3 THEN Ya= Ya/FNHsmooth(Xa. ltd) 
X=LGT(Xa/Xmin)-Sfx 
Y=LGT< Ya/Ynun)*Sfy 
IF Sym>3 THEN PRINT ”SM” 

IF Sym<4 THEN PRINT “’SR 1,4, 2. 4” 

PRINT ”PA”.X,Y,”” 

IF Sym>3 THEN PRINT "SR 1.2, 1.6” 

IF Sym=4 THEN PRINT ”UC2 . 4 . 99 . 0 , -8 . -4 , 0 , 0 . 8 , ^ . 0 : " 

IF Sym = 5 THEN PRINT "UC3 , 0 , 99 , -3 , -6 . *3 . 6 , 3 , 6 . 3 . -6 : " 
IF Sym=6 THEN PRINT "UCO . 5 , 3 . 99 . 3 . -8 . -6 , 0 . 3 , 8 : " 

IF Sym = 7 THEN PRINT "UCO , -5 . 3 . 99 . -3 , 3 . 5 . 0 . - 3 . -8 : " 

NEXT I 

BEEP 

ASSIGN aiFile TO - 
GOTO 5690 
END IF 

PRINT "PU SM" 

BEEP 

INPUT "WANT TO PLOT A POLYNOMIAL ( I = Y . 0=N ) ?" . Go_on 

IF Go_on=1 THEN 

BEEP 

INPUT "ENTER LOWER AND UPPER X-LIMITS" , X 1 1 , X ia 
FOR Xx=0 TO Cx STEP Cx/200 
Xa=Xm i n* 1 0 ' Xx 

IF Xa<XiI OR Xa>Xla THEN 6380 

Y a=FNPo iy ( Xa) 

Y = LGT< Ya/Ynun)^Sfy 
X=LGT(Xa/Xmin)-Sfx 
IF Y<0 THEN Y=0 
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D JbU 

6370 

8380 

6390 

8400 

6410 

8420 

6430 

6440 

6450 

6460 

6470 

6480 

6490 

6500 

6510 

6515 

6520 

6525 

6530 

6535 

6540 

6545 

6550 

6555 

6560 

6565 

6570 

6575 

6580 



ir Y>IUU "HEN iiDTO 6o8U 
PRINT ••PA" .X. Y .“PD“ 

NEXT <x 
END 

PRINT "PU PA 0.0 SPO" 

SUBEMD 

OEF FNHsmoo th ( X . I td ) 

Hs-FNPoLy(X)/X 

IF Itd=l THEN Hs=Hs*. 83347 

RETURN Hs 

FNEND 

DEF FNPoly(X) 

Poiy--4.4l23718E+2-X*(6,8123917E^2-X*3.7416863E'^2) 

RETURN Poly 
FNEND 

DEF FNPvst< Tsteam) 

DIM K<8) 

DATA -7.691234564.-26.08023696.-168. 1706546,64.23285504,-1 18.9646225 
DATA 4.16711732,20.9750676.1 .E9.S 
READ K(*) 

T = <Tsteain+273. 15)7647.3 
Sujn = 0 

FOR N=0 TO 4 

Sum=Sam+K(N)»( l-T) <N+1 > 

NEXT M 

8r = Sojn/<T*< l'^K<5)-< l-T) ^-K<6)^< l-T) 2) )-(1-T)/(K<7)-<l-T> 2 ^K<8>) 
Pr-EXP(Br) 

P=2212000U*Pr 
RETURN P 
FNEND 
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APPENDIX B 



SAMPLE CALCULATIONS 



Data set number fifteen of file GT6H29 corresponding to 
tube number 6 was chosen for the sample calculation. 



A. 



TEST-SECTION 


DIMENSIONS 




O.D. 


= 


25.30 


(mm) 




= 


1.05 


(mm) 


fins/mm 


= 


1.02 




inter-fin gap 


= 


0.25 


(mm) 


°2 


= 


23.10 


(mm) 


°1 


= 


17.86 


(mm) 


D 

o 


= 


22.22 


(mm) 


D. 

1 


= 


20.32 


(mm) 


^u 


- 


31.75 


(mm) 


L 


= 


50 . 80 


(mm) 



B. MEASURED PARAMETERS 



RES 



= 70.0 (volts) 

= 4.8 (volts ) 

= 2.031 (ohms) 

= 383 (W/m.K) 

= 54.23 (°C) 

= 54.29 (°C) 

= 53.91 (°C) 
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53.65 (°C) 




= 53.40 (°C) 

= 53.31 (°C) 

= 53.22 (°C) 

= 53.54 (°C) 

= 47.51 i°C) 



PROPERTIES OF FREON- 11 3 AT SATURATION TEMPERATURE 

,-4 



y 

P 

C. 



= 4.948 X 10 (kg/m.s) 

= 1510.0 (kg/m^) 

= 0.973 (kJ/kg.K) 

= 0.0715 (W/m.K) 



D. HEAT-FLUX CALCULATION 

The heat transfer rate from the cartridge heater: 



= VI (W) and I = V^^^g/R (amps) 

I = 4.8/2.031 = 2.363 (amps) 



Q„ = 70.0 x 2.363 = 165.410 (WO 

H 

The cross-sectional area of machined ends: 

A = tt/ 4.0(D^ - D^) (m^) 

^ 0 3 . 

A^ = tt/ 4.0 (0.02222^ - 0.02032^) = 63.481 xio"^ (m^) 
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The tube outside wall perimeter: 



P = ttD = 3.1416 x 0.02222 = 6.981 x 10 ^ (m) 

0 

The corrected length of ends: 

Lc = + t/2.0 t - - D. 

L = 0.03175 + (0.02222 - 0.02032)/2.0 = 32.7 x 10'^ (in) 

Other thermophysical properties: 

V = y/p = 4.948 X 10“^/1510.0 = 3.277 xio""^ (m^/s) 

a = k/pCp = 0.0715/(1510.0 x 973) = 4.866 xio"^ (m^/s) 

1 Ap _ 1 ,1509.77 - 1510.26, 

^ p AT 1510^ 47.61 - 47.41 ^ 

= 0.00162 (1/K) 

Pr = v/a = 3.277 X 10"’^/4. 866 X 10"^ = 6.734 



Temperature at the base of the T-shaped fins: 



T 



B 



ln(D /D ) 

T - o = — 

AVG 2 7T L 



T 



T 



AVG 



B 



8 

= ( I T^)/8.0 

n=l 

53.69 - 165.41 



= 53.69 (°C) 



ln(23.1/17.86) 
2 ttx 0. 0508x383 



53.34 (°C) 
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53 . 342 



47.51 



5.832 (K) 



B 



= Tg-T 



'SAT 



Average heat-transfer coefficient at the unenhanced 
ends : 



h 



J5_ 

D 

o 



0.60 +0.387 



r 3 f hP 1/2 ,, 


1/6 


, hP , 1/2 

A ^ J 

C C 




[1 + (0.559/Pr)^/l^]®/2'^ 



\ 

h = 187.49 (W/m^.K) 



/ 



Heat- transfer rate through unenhanced ends: 

Qp = vh P A^ 0g Tanh mL^ 

, hP ,1/2 ,187.489 x 6.981 x lo"^ , ^ ,„3 

m = ( . V ■) = ( = 1.698 X 10 

C C 0.0715 x 63.481 x 10 ° 

Qp = 3.29 (W) 

QlOSS = 



Heat flux through Gewa-T surface: 



Q = Qh - Q 



•LOSS 
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Q 



165.41 



6.579 



158.83 (W) 



q = Q/Ag 

Ag = TT X 23 . 1 X lo"^ X 50 . 8 X 10~^ = 3.687 xio"^ (m^) 



q = 



158.83 



3.687 X 10 



-3 



43.978 X 10^ (W/m^) 



The following are the results obtained from the computer 
by running the data reduction program (DRP) (see Appendix A) 

Q = 44.466 10 (W/m^) 

0g = 5.76 (K) 

h- = 187.39 (W/m^.K) 

The small difference in results obey the inaccuracy of 
hand calculations with the use of about three significant 
digits only. 
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APPENDIX C 



UNCERTAINTY ANALYSIS 



The same data set (set No. 15 of file GT6H29) that was 
used for the sample calculation was chosen for the uncertainty 
analysis; therefore dimensions of test section, measured and 
calculated parameters found in sample calculation were used 
in this analysis. All uncertainties are presented as a 
percentage of the calculated parameter. 

A. UNCERTAINTY IN SOURCE HEAT-TRANSFER RATE 



Q 



VI (W) 



I 



Vres/^^ (amp) 



6V 



6V 



RES 



0.05 (volts) 



where : 



6 = uncertainty in measurement and calculation 



61 



6V„^^/R = 0.05/2.031 = 0.025 (amp) 

KijO 





V I 





1.06 percent 
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B. UNCERTAINTY IN SURFACE AREA 



B 



u D2L 



6D2 = 5 ^ 



= 0.1 






0.1 2 0.1 2 

^^ 23 . 1 ^ ^ 50 . 8 ^ ^ 



= 0 



C. UNCERTAINTY IN AT 



AT = T - T 6T = 6T 

^ W SAT W ° SAT 



6AT 

AT 



2 5T 2 1/2 

[ ( + (- ] 



AT 



AT 



5AT 

AT 






1/2 



D. UNCERTAINTY IN HEAT FLUX 



q = 



Qh - 2Qp 






. 5Q 2 25Q 2 ^^1 

M = [( H 2 ^2 ] 

q -2Q ^ ^ ^ ^ A. 






In this case, Q = 165.41 (W) and Q. 

ll 



(mm) 



48 percent 



0.15 (°C) 



3.43 percent 



2 1/2 
) ] 



= 3.29 (W) 
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Assuming the same 



Therefore, Q = 2.0 percent of Q„. 

r n 

proportion in the uncertainty for Q : 

F 



6Q„ = 0.02 6Q„ = 0.02x1.753 = 0.035 (W) 

r rl 



^ ^ ^ + (-0.0048)^] 



1/2 



158.83 



^ = 



1.26 percent 



E. UNCERTAINTY IN BOILING HEAT-TRANSFER COEFFICIENT 



h = ^ 



AT 



6h 

h 



[(M)2 + (. «AT) 

L V V ^ AT ^ 



1/2 



— ]_/2 
— = [(0.0126)^ + (-0.0343)^] = 3.65 percent 
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